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The next core-collapse supernova in our Galaxy will be a spectacular event, with 
some 10"* neutrino detections in total expected among several detectors. This 
data will allow unprecedented tests of neutrino properties and new opportunities 
in astrophysics. In this paper, I focus on two main topics: (1) Measurement of 
the Ufj, and i/t masses by time-of-flight, with an emphasis on introducing as little 
supernova model dependence as possible, and (2) Methods for locating a supernova 
by its neutrinos in advance of the light, which may allow improved astronomical 
observations. In the latter, I also discuss the recent result that the positrons from 
-h p ^ e'^ n are not isotropically emitted, as commonly thought. 



1 Supernova Neutrino Mass Measurements 

1.1 Current Knowledge About Neutrino Masses 

What is known about the neutrino masses, in terms of direct experimental 
evidence? After several decades of experiments on the tritium beta spectrum, 
the miblished limit on the mass of the (and hence also i>e) is now about 
3 eVEl. These results are still attended by some controversy, due to apparent 
systematic errors near the endpoint (discussed in detail at a recent workshopH 
at the Institute for Nuclear Theory in Seattle). Much less is known about the ia,, 
mass (< 170 keV from tt decay) B, and the Vr mass (< 18 MeV from r decay) □. 
It would be a great advance if the latter two limits could be improved to the 
eV or tens of eV level, and it seems unlikely that any terrestrial experiment 
could do that. 

There are two indirect arguments which may restrict the z/^ and Vt masses. 
The first says that the predicted background of relic neutrinos (« 100/cm'^) 
should not overdose the universe. This leads to a limit on the sum of the 
neutrino masses of about 100 eV in the most conservative case, and about 
10 eV if reasonable valuesJor the non-baryonic dark matter density and the 
Hubble constant are usedlj. The second is based on the recent evidence for 
neutrino oscillations (sensitive only to differences of neutrino masses). For 
three neutrino flavors, there are two independent mass differences, which can 
be chosen to explain the solar and atmospheric neutrino data. The overall 
scale can be fixed by using the limit from tritium beta decay, and since the 
observed mass differences are small, then all masses are below a few eVEl. 
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These indirect arguments sound compelling, but they are not hard to 
evade. The cosmological bound doesn't apply if the neutrinos are allowed 
to decay. There^ are many such models for heavy v-t masses, motivated by 
particle physics Q or astrophysics eI. The oscillation bound doesn't apply if 
there are more free mass differences than there are positive signals of neutrino 
oscillations. For example, suppose the atmospheric neutrinos are undergoing 
maximal ^ Vr mixing, LSND is ruled out, and the solar neutrinos are un- 
dergoing <-> Vsteriie mixing. Then nothing prevents the 1^2 and masses 
from having a tiny splitting near 50 or even 500 eV, and leaving the mass 
below a few eV. 

It may be possible in the future to make newj-astrophysical tests of the 
neutrino masses, using data on. large-scale structurela't2l, the cosmic microwave 
background til, weak lensinglla, and the Lyman-alpha forest 113. While these 
techniques claim sensitivities of order 1-10 eV for the sum of the neutrino 
masses, none are based on direct detection of neutrinos, and hence may be 
vulnerable to uncertainties in their assumptions. 

Thus while indirect evidence will be valuable, there is still a need for a 
direct measurement of the and masses (or, in the presence of mixing, the 
masses of the heavy mass eigenstates). 

1.2 Supernova Neutrino Emission 

When the core of a large star (M > SMq) runs out of nuclear fuel, it collapses, 
with a change in the gravitational binding energy of about 3 x 10^^ ergs. This 
huge amount of energy must be removed from the proto-neutron star, but the 
high density prevents the escape of any radiation. Inside of the hot proto- 
neutron star, neutrino-antineutrino pairs of all flavors are produced. Despite 
their weak interactions, even the neutrinos are trapped and diffuse out over 
several seconds, in the end carrying away about 99% of the supernova energy. 
(I have neglected the flux of Ve neutrinos necessary to change the core from 
N Z nuclei into a neutron star, because the total energy release of this flux 
is of order 1% of the pair emission phase). 

When the neutrinos are about one mean free path from the edge, they 
escape freely, with a thermal spectrum (approximately Fermi-Dirac) charac- 
teristic of the surface of last scattering. Because different flavors have different 
interactions with the matter, and because the neutron star temperature is 
decreasing with increasing radius, the neutrino decoupling temperatures are 
different. The i/^ and Vr neutrinos and their antiparticles have a temperature 
of about 8 MeV, the De neutrinos about 5 MeV, and the I'e neutrinos about 3.5 
MeV. Equivalently, the average energies are about {E) ~ 25 MeV, ~ 16 MeV, 
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and ~ 11 MeV, respectively. The luminosities of the different neutrino flavors 
are approximately equal at all times. The neutrino luminosity rises quickly 
over 0.1 s or less, and then falls off over several seconds. The SN1987A data 
can be reasonably fit to a decaying exponential with time constant r = 3 s. 
The detailed form of the neutrino luminosity used below is less important than 
the general shape features and their characteristic durations. 

The estimated core-collapse supernova rate in the Galaxy is about 3 times 
per century ta. The present neutrino detectors can easily observe a supernova 
anywhere in the Galaxy or its immediate companions (e.g., the Magellanic 
Clouds). Unfortunately, the present detectors do not have large enough vol- 
umes to observe a supernova in even the nearest galaxy (Andromeda, about 
700 kpc away). 

1.3 Time- of -Flight Concept 

Even a tiny neutrino mass will make the velocity slightly less than for a massless 
neutrino, and over the large distance to a supernova will cause a measurable 
delay in the arrival time. A neutrino with a mass m (in eV) and energy E (in 
MeV) will experience an energy-dependent delay (in s) relative to a massless 
neutrino in traveling over a distance D (in 10 kpc, approximately the distance 
to the Galactic center) of 



where only the lowest order in the small mass has been kept. A distance of 10 
kpc corresponds to a travel time of about 10^^ s, and the delays of interest are 
less than 1 s. 

If the neutrino mass is nonzero, lower-energy neutrinos will arrive later, 
leading to a correlation between neutrino energy and arrival time. This is 
exploitable in some of the charged-current detection reactions for and Dg, 
since the neutrino energy can be reasonably measured. Using this idea, the 
next supernova will allow sensitivity to a Ve iipass down to about 3 eV 113. 
However, the terrestrial experiments now claimEl to rule out a mass that large, 
so I ignore the Ve mass in describing the supernova neutrino signal. 

For measuring the masses of the v^^ and v^- neutrinos, one must realize 
that only neutral-current detection reactions are possible (the charged-current 
thresholds are too high) . Since the neutrino energy is not measured in neutral- 
current interactions, an energy-time correlation technique cannot be used (the 
incoming neutrino energy is not determined since a complete kinematic recon- 
struction of the reaction products is not possible). 




(1) 
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Instead, the strategy for measuring the Vt mass is to look at the difference 
in time-of-flight between the neutral-current events (mostly u^^Vt,D^^ and 
because of the higher temperature of those flavors) and the charged-current 
events (just and Ve)- We assume that the is massless and will ask what 
limit can be placed on the v-r mass (other cases will be discussed below). 
There are three major complications to a simple application of Eq. (|l|): (i) 
The neutrino energies are not fixed, but are characterized by spectra; (ii) 
The neutrino pulse has a long intrinsic duration of about 10 s, as observed 
for SN1987A; and (iii) The statistics are finite. In some early worktEl, the 
duration of the pulse was thought to be of order 100 milliseconds, nearly a 
delta-function. Here (considering the case of the smallest detectable mass), 
the delay is much less than the width of the pulse, which is what makfis.ptiie 
problem much more difficult (the large-mass case is covered elsewhere Il3il3). 
Given that the energy-time correlation method cannot be used, we must turn 
to considerations of the shape of the event rate as a function of time. 



1.^ Neutrino Scattering Rate 



For a time-independent spectrum f{E), and luminosity L{ti), the double- 
differential number distribution at the source is 



dE^dti 



(2) 



The neutrino spectrum and average energy (E^) are different for Cjoch flavor, 
but the luminosities L{ti) of the different flavors are roughly equal IIj, so that 
each flavor carries away about 1/6 of the total binding energy Eb- For a 
massive neutrino, the double-differential number distribution at the detector 
is 

d'^N^ L{t- At{Ey)) 



dE^dt {E^) 
and the scattering rate for a given reaction is 



dN^ 
dt 



C / dE,f{E, 



a{E,) 



L{t-At{E,)) 
Eb/6 



The overall constant for an II2O target is 



C = 9.21 



Eb 



ergs 



1 MeV 



T 



10 kpc 



D 



det. mass 



1 kton 



(3) 



(4) 



(5) 



where n is the number of targets per molecule for the given reaction. For a 
D2O target, the leading factor becomes 8.28 instead of 9.21. 
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Table 1: Calculated numbers of events expected in SK (32 kton H2O) with a 5 MeV threshold 
and a supernova at 10 kpc. The other parameters (e.g., neutrino spectrum temperatures) 
are given in the text. In rows with two reactions listed, the number of events is the total for 
both. The second row is a subset of ths first row that is an irreducible background to the 
reactions in the third and fourth rowscj. The detected final-state particles are e+,e~,7. 



reaction 


number of events 


De + P ^ + n 


8300 


De+P^e+ + n {E^+ < 10 McV) 


530 




355 


Z/^ +1« ^ t/r + 7 + ^ 

+16 z?^ + 7 + X 


355 


Ve + e~ ^ Ue + e~ 
De + e~ Ve + 


200 


+ ^ + 
t'/i + ^ + 


60 


U-r + e~ ^ Vr + eT 
Vr + e~ Vr + 


60 



For a massless neutrino, At{E,^) = 0, so the luminosity comes outside 
of the integral as L{t) and completely specifies the time dependence. For a 
massive neutrino, the time dependence from the luminosity is modified in an 
energy-dependent way by the mass effects, as above. 

1.5 Separation of the Massive Signal 

The scattering rates of the massless neutrinos {v^, or v^) can then be used to 
measure the shape of the luminosity as a function of time. The scattering rate 
for the massive neutrinos (vr) can be tested for additional time dependence 
due to a mass. We define two rates: a Reference R{t) containing only massless 
events, and a Signal S{t) containing some fraction of massive events. 

The Reference R{t) can be formed in various ways, for example from the 
charged-current reaction Ve + p ^ + n in the light water of SK or SNO. 
(The numbers of events expected for all reactions are given in Table I for SK 
and Table II for SNO). nnn 

The Signal S{t) can be based on various nfiutral-current reactions EZlEflO, 
though here I will focus on the results for SNOES. Thus the primary component 
of the Signal S{t) are the 485 neutral-current events on deuterons. With the 
hierarchy of temperatures assumed here, these events are 18% {ve + Ve), 41% 
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Table 2: Calculated numbers of events expected in SNO for a supernova at 10 kpc. The 
other parameters (e.g., neutrino spectrum temperatures) are given in the text. In rows with 
two reactions listed, the number of events is the total for both. The middle column indicates 
the detected particle(s). The notation v indicates the sum of v^, i^^, and v^, though tiiey 
do not contribute equally to a given reaction, and X indicates either n+'^^O or p+^^N. 



Events in 1 kton D2O 




n 


485 




n 




i^e + rf^e~+p + p 


e~ 


160 


i^e + cJ^e^ + n + n 


e^nn 




z/+^'*0^!/ + 7 + X 


7, 7n 


20 


p+i6o^p + 7 + X 


7, 771 






n 


15 


j> +16 ^ z/ + n 


n 








10 


V + er v -\- 


er 




EvtHits in 1.4 Ivton II2O 


Ve + P ^ e'^ + n 


e+ 


365 




1 


30 


p+16 0^zy + 7 + X 


7 




V + e" v + e~ 


e~ 


15 


p + e~ D + e~ 


e~ 





{Vfji+Vfj), and 41% {vt-\-Vt). The flavors of the neutral-current events of course 
cannot be distinguished. Under our assumption that only Vt is massive, there 
is aheady some unavoidable dilution of S{t). 

In Fig. 1, S{t) is shown under different assumptions about the u^- mass. 
The shape of R{t) is exactly that of S{t) when m^^ = 0, though the number 
of events in R{t) will be different. The rates R{t) and S{t) will be measured 
with finite statistics, so it is possible for statistical fluctuations to obscure the 
effects of a mass when there is one, or to fake the effects when there is not. We 
determine the mass sensitivity in the presence of the statistical fluctuations 
by Monte Carlo modeling. We use the Monte Carlo to generate representative 
statistical instances of the theoretical forms of R{t) and S{t), so that each run 
represents one supernova as seen in SNO. The best model-independent test 
of a Ur mass seems to be a test of the average arrival time {t). Any massive 
component in S{t) will always increase {t), up to statistical fluctuations. 
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Figure 1: The expected event rate for the Signal S{t) at SNO in the absence of fluctuations 
for different i/t masses, as follows: solid line, eV; dashed lines, in order of decreasing height: 
20, 40, 60, 80, 100 eV. Of 535 total events, 100 are massless {ue + Pg), 217.5 are massless 
(Ufi + Pp), and 217.5 are massive {ur + Pt)- These totals count events at all times; in the 
figure, only those with t < 9 s are shown. 
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1.6 {t) Analysis 

Given the Reference R{t) (i.e., the charged-current events), the average arrival 
time is defined as 

{t)H = ^ , (6) 

where the sum is over events in the Reference. The effect of the finite number 
of counts Nr in R{t) is to give {t)R a statistical error: 

^mn) = :^^^M. (7) 

For a purely exponential luminosity, {t)ji = \J {t'^)R — = t ~ 3 s. Given 
the Signal S{t) (i.e., the neutral-current events), the average arrival time {t)s 
and its error 5{{t)s) are defined similarly. The signal of a mass is that the 
measured value of {t)s — {t)R is greater than zero with statistical significance. 

Using the Monte Carlo, we analyzed 10^ simulated supernova data sets for 
a range of Vt masses. For each data set, {t)s — {t)R was calculated and its 
value histogrammed. These histograms are shown in the upper panel of Fig. 2 
for a few representative masses. (Note that the number of Monte Carlo runs 
only affects how smoothly these histograms are filled out, and not their width 
or placement.) These distributions are characterized by their central point and 
their width, using the 10%, 50%, and 90% confidence levels. That is, for each 
mass we determined the values of {t)s — {t)R such that a given percentage of 
the Monte Carlo runs yielded a value of {t)s — {t)R, less than that value. With 
these three numbers, we can characterize the results of complete runs with 
many masses much more compactly, as shown in the lower panel of Fig. 2. 
Given an experimentally determined value of {t)s — (i)_R, one can read off the 
range of masses that would have been likely (at these confidence levels) to have 
given such a value of {t)s — (i)_R in one experiment. Prom the lower panel of 
Fig. 2, we sec that SNO is sensitive to a I'r mass down to about 30 eV if the 
SK R{t) is used, and down to about 35 eV if the SNO R{t) is used. 

We also investigated the dispersion of the event rate in time as a measure 
of the mass. A mass alone causes a delay, but a mass and an energy spectrum 
also causes dispersion (as does the separation of the massive and massless por- 
tions of the Signal). We defined the dispersion as the change in the width 
\/ {t'^)s ~ ~ \/ {i'^)R — {t)^R- We found that the dispersion was not statisti- 
cally significant until the mass was of order 80 eV or so, i.e., a plot very similar 
to Fig. 2 can be formed for the width, and the error band is much wider. Since 
the dispersion is irrelevant, the average delay is well-characterized by a single 
energy, which for SNO is Ec ~ 32 MeV. 



8 



-I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — l- 



m = OeV 



m = 50 eV 



m = 100 eV 




-0.2-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



<t>s-<t>R 



[s] 



Figure 2: The results of the (t) analysis for a massive Vt, using the Signal S{t) from SNO 
defined in the text. In the upper panel, the relative frequencies of various (t)s — {'t)R values 
are shown for a few example masses. The solid line is for the results using the SK Reference 
ii(t), and the dotted line for the results using the SNO R(i)- In the lower panel, the range of 
masses corresponding to a given {t)s — {t)R is shown. The dashed line is the 50% confidence 
level. The upper and lower solid lines are the 10% and 90% confidence levels, respectively, 
for the results with the SK R(t). The dotted lines are the same for the results with the SNO 
R{t). 
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1.7 Analytic Results 

Another nice feature of the proposed {t) analysis, besides its simpUcity, is that 
good estimates can be made analytically. The characteristic delay is 

(t) s - (t)i?^frac(m>0)x 0.515 f—j D, (8) 

where frac(m > 0) is the fraction (about 1/2) of massive events in the neutral- 
current Signal. The above numerical analysis shows that the mass effects 
occur most significantly in the delay and not the width (are characterized 
with a single energy, and negligible dispersion) . That characteristic energy Ec 
is roughly at the peak of J{E)a{E) (see Fig. 3). This formula for the first 
moment is always true; the point of the Monte Carlo analysis was to show that 
the other moments are not changed, so this completely describes the data. 

If the cross section (j{E^) depends on energy as E" (a ~ 2 for + d), 
then the characteristic energy Ec ^ {2 + a)T and the thermally-averaged cross 
section is proportional to T", where T is the t'^ and i>t temperature. For 
a = 2, and a neutron detection efficiency e„, the following scaling relations 
hold: 

{t)s-{t}B.-{j;) D. (9) 




For a non-exponential luminosity, t is more generally the width of the pulse. 
If zero delay is measured, then the mass limit is determined by the largest 
positive delay that could have fluctuated to zero, and thus 

™H™*-T3/4rV2,-i/4. (12) 

Note that is independent of D for any supernova distance in the Galaxy (i.e., 
as long as there is a reasonable number of counts). Note that the mass limit 
scales with the fourth root of the number of events (via the efficiency e„). 

Besides the scaling, Jjlie value of mum at fixed values of these quantities 
can also be found easily E3. The anahdjic results above are in excellent agree- 
ment with the full numerical results 113. Over reasonable ranges of the input 
parameters, the results hardly change— However, some estimates of the v^- mass 
sensitivity for proposed detectors E3'a assumed a very short duration of the 
supernova (r ~ 0.5 s). If such a short duration were valid, then the SNO 
sensitivity would also be of order 10 eV. 
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Figure 3: The shapes of response functions for different reactions. The sohd Une is a Fermi- 
Dirac distribution /(Ei/) with temperature 8 MeV. The dashed hnes are products of that 
with the various neutral-current cross sections, i.e., f{Ev)a{Ei,). The normahzations were 
chosen so that all curves would have the same height. From left to right, the cross sections 
are Jf+e", v + d, v+^^C, i^+plO. For ^^c, only the level at 15.11 MeV was used. For ^^O, 
a fitliJ to the full calculation CJ was used. In a given experiment, only the integral of one 
of these curves is measured, but each samples the underlying f{Eu) differently. The peak 
positions are what is called Ec in the text. 



1.8 How Model- Dependent are the Results? 

Above, I discussed the proposed {t) technique for analyzing the supernova 
neutrino data for a Ut mass. This technique was designed to be as model- 
independent as possible. One underlying, theme is to calibrate the assumed 
models as much as possible from the datatZHla. 

While is is probably a good approximation that the luminosities of the 
different flavors are equal as a function of time, the normalization of the lumi- 
nosity drops out in the definition of (i), and hence the results depend on the 
weaker assumption that the shapes of the luminosities are approximately equal 
as a function of time. The differences among the different flavors at the very 
early stages (for example, the higher Ve flux) persist over a duration much less 
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than 10 s, and hence have no effect on the analysis. 

Another very important point is that the details of the shape of the lu- 
minosity as a function of time are never used, because the delay is measured 
just from the {t) of the neutral-current events minus the same for the charged- 
current events. Why is this enough? In general, one could describe the scatter- 
ing rates of charged-current and neutral-current events by their normalizations 
and a series of all moments {{t), {t^), etc.). While this would likely require 
many terms, it is efficient to consider the difference in the moments between 
the two. The mass effects do not change the normalizations (the numbers of 
events), and we have shown that the higher moments (e.g., the width) are 
not significantly changed by a mass. That is, the only statistically significant 
effects of a mass on the event rate difference are in the average arrival time 
(<) . If the details of the supernova model are not known from theory, then {t) 
is what is called a sufficient statistic for the masscJ. 

The only aspect of the shape of the luminosity L{t) that was used in the 
analysis was its duration, and this appears only in the error on the delay. 
Thus the particular assumption of using an exponential luminosity completely 
disappears from the problem. Any other reasonable model of similar duration 
would give a very similar result. 

One also needs to know the temperature to determine the mass limit. 
Since none of the neutral-current reactions measure the incoming neutrino en- 
ergy, there is no direct spectral information. The Vt temperature must be 
determined from the number of events. Since different neutral-current reac- 
tions have different energy dependence, one can crudely reconstruct the un- 
derlying spectrum (see Fig. 3). The moments analysis of the rates reveals that 
only the average delay is significant, and thus that just a single energy (the 
characteristic energy Ec) is relevant, instead of the fine details of the spectral 
shape. 

The final results for the mass sensitivity are given in Table III. 
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2 Supernova Location by Neutrinos 

2. 1 Overview 

There has been great interest recently in the question of whether or not a 
supernova can be located by its neutrinos. If so, this may offer an opportunity 
to give an early warning to the astronomical community, so that the supernova 
light curves can be observed from the earliest possible time. An international 
supernova early alert network has been formed for this purpose, and the details 
of its implementationE^uj were the subject of a recent workshopEZl. One of the 
primary motivations for such a network is to greatly reduce the false signal rate 
by demanding a coincidence between several different detectors. The second 
motivation is to locate the supernova by its neutrino signal. 

The interest in the latter is driven by two considerations. First, the neutri- 
nos leave the core at the time of collapse, while the electromagnetic radiation 
leaves the envelope some hours later, depending on the stellar mass. Thus the 
neutrino observations can give the astronomers a head start. Second, the next 
Galactic supernova will likely be in the disk, and hence obscured by dust, and 
perhaps never visible optically. 

There are two classes of techniques to locate a supernova by its neutrinos. 
The first class of techniques is based on angular distributions of the neutrino 
reaction products, which can be correlated with the neutrino direction. In 
this case, a single experiment can independently announce a direction and its 
error. The second method of supernova location, triangulation, is based on the 
arrival-time differences of the supernova pulse at two or more widely-separated 
detectors. This technique would require significant and immediate data sharing 
among the different experiments. 

2.2 Angular Distributions 

Neutrino-Electron Scattering: Electron Angular Distribution 

In neutrino-electron scattering, the scattered electrons have a very forward 
angular distribution, due to the small electron mass. At supernova neutrino 
energies, the angle between the incoming neutrino and the outgoing electron 
is about 10°, depending somewhat on neutrino energy and flavor. However, 
multiple scattering of the struck electron smears out its Cerenkov cone, washing 
out the dependence on energy and flavor, and introducing an angular resolution 
of about 25°. 

Naively, if the one-sigma width of the electron angular distribution is 25°, 
then the precision with which its center (i.e., the average) can be defined given 
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Ns events is 




where Ns is the number of events. For SK EH, Hs-, — 320, so the cone center 
could be defined to within about 1.5°. For SNOEfl (using both the hght and 
heavy water), Ns — 25, so the cone center could be defined to about 5°. The 
equivalent error on the cosine is 6{cos6) ~ {69)^/2, i.e., 3 x 10"'' and 4 x 10"'^, 
respectively. 

However, the neutrino-electron scattering events are indistinguishable on 
an event-by-event basis from other reactions with only a single electron or 
positron detected. The primary background in SK or the light water in SNO 
is thus Pe+P — e^"+n (in the heavy water in SNO, there are several background 
reactions) . 

This problem is similar to the SK solar neutrino studies, in which the 
neutrino-electron scattering events have to be separated from the intrinsic de- 
tector background (due to radioactivity, etc). The SK results are presented as 
an angular distribution in the cosine along the known direction to the Sun, and 
clear bump in the solar direction is seen, with signal/noise ~ 1/10. Though the 
concept is the same, for the supernova events the signal/noise is only ~ 1/30. 
Taking this into account c3 makes the centroiding error larger by a correction 
factor of ~ 4 for both SK and SNO. With some cuts on energy, it may be 
possible to reduce this correction factor to ~ 2 — 3. 

Neutrino-Nucleus Scattering: Lepton Angular Distributions 

In the charged-current reactions of I'e and Vg on nucleons or nuclei, the outgoing 
electrons or positrons typically have angular distributions of the form 

~ l-HWea(£;^)cos6l, (14) 

dcost^ 

where 9 is the angle between the neutrino and electron (or antineutrino and 
positron) directions in the laboratory (where the nuclear target is assumed to 
be at rest) and Ve is the lepton velocity in c = 1 units. At energies higher 
than for supernova neutrinos, terms proportional to higher powers of cos9 
appear, and at the highest energies, all reaction products are strongly forward 
simply by kinematics. Just above threshold, where the lepton velocity is small, 
the angular distribution obviously becomes isotropic. At supernova neutrino 
energies, = 1. 

Typically, the asymmetry coefficient a{E^) is not large, so the angular 
distribution is weak, but the number of events may be large (for i>e+P e^+n, 
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there are nearly 10^ events expected in SK). Thus it may be possible to use 
these events to locate the supernova. We first consider how well one could 
localize the supernova assuming that a is known and constant. 

Given a sample of events, one can attempt to find the axis defined by 
the neutrino direction. Along this axis, the distribution should be fiat in the 
azimuthal angle (f) and should have the form in Eq. in cos6'. Along any 
other axis, the distribution will be a complicated function of both the altitude 
and azimuthal angles. We assume that the axis has been found numerically, 
and ask how well the statistics allow the axis to be defined. A convenient way 
to assess that is to define the forward-backward asymmetry as 

where Np and A^^ are the numbers of events in the forward and backward 
hemispheres. The total number of events is A^ = Np + Nb- Note that ApB 
will assume aiLextremal value A"^^ along the correct neutrino direction. It 
can be showncj that 



SApB 

/7VV V2/ 

where the error is nearly independent of a for small |a|, which is the case for 
the reactions under consideration. 

In the above, the coordinate system axis was considered to be correctly 
aligned with the neutrino direction. Now consider what would happen if the 
coordinate system were misaligned. While in general, all three Euler angles 
would be needed to specify an arbitrary change in the coordinate system, 
symmetry considerations dictate that the computed value of ApB depends 
only upon one - the angle 9 between the true and the supposed neutrino axis. 
Thus ApB is some function of 9 if the axis is misaligned. Using a Legendre 
expansion, one can show that 

ApB{e) ^^cos0. (17) 
The error on the alignment is then 

2 2 1 

S{cos9) ^ —rSApB Tl^- (18) 
\a\ \a\ y/N 

Treating the nucleons as infinitely heavy, the coefficient a in Eq. ([l^ ) is 
related to the competition of the Fermi (no spin flip) and Gamow- Teller (spin 
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flip) parts of the matrix element squared: 

_ \Mf\'-\Mgt? 

|M;^|2 + 3|A/gtP ■ ^ ^ 

Naive use of this formula for S'e+P — > e+ + n gives a = —0.10 for {\Mgt/Mp\ = 
1.26). One expects 10^ events in SK, so that d{cos9) ~ 0.2. Even though the 
asymmetry parameter a is quite small, the number of events is large enough 
that this technique in SK could give a reasonable pointing error (SNO only 
has 400 of these events, so S{cosd) ~ 1.0, which is too large to be useful). 
For the charged-current reactions on deuterons, De + d ^ e+ + n + n and 
i/e + d ^ + p + p, one would have a = —1/3 {Mp /Mgt — 0). If these 
channels could be combined (160 events in total), then 5{cos9) ~ 0.5, which is 
again rather large. 

In general, the coefficient a(E^) has energy dependence coming from recoil 
and weak magnetism corrections, each of order 1/M, jiihexe.M is the nucleon 
mass (such terms also change the total cross section Ej't3Eil, which has some 
important applications e3). For the reaction Pf.+p e+ +ri, these corrections 
have a dramatic effect on the angular distribution, making it backward at low 
energies, isotropic at about 15 MeV, and forward at higher energies. Averaged 
over the expected 9^ spectrum from a supernova, one obtains (cos0) = +0.08, 
i.e., a forward distribution. Using the common assumption that the angular 
distribution is backward (with a = —0.10) would lead one to infer exactly the 
wrong direction to the supernova. 

It is convenient to describe the angular distribution by the average cosine, 
weighted by the differential cross section. In the limit that Eq. ( p^ holds, 

{cos0) = ^Vea{E.). (20) 

For + p — > e"*" + n, the variation of the average CQsine can be written as 
(keeping only the largest terms of the full expressionEil) : 

(cos6') ~ -0.034we + 2.4-^. (21) 

The factor 2.4 is actually 1.0 + 1.4, the first term from recoil, which always 
makes the angular distribution more forward, and the second from weak mag- 
netism (which would change sign for the reaction v,, + n ^ e~ + p, if there 
were free neutron targets). The variation of {cos 6) with E^, is shown in Fig. 4. 

Since the deuteron is weakly bound, the recnil and weak magnetism cor- 
rections to (cos 6) can be reasonably estimated tiJ. The effects of recoil and 
weak magnetism add in the i>e + d channel, and partially cancel in the + d 
channel, as seen in Fig. 5. 
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Figure 4: Upper panel: total cross section for t'e + p — > e"*" + n; bottom panel: (cos 9) for the 
same reaction; both as a function of the antineutrino energy. The solid line is our 0(1/M) 
result and the short-dashed line is the 0(1) result. The long-dashed line is the high-energy 
formula of Llewellyn-Smith, and the dot-dashed line contains our threshold modifications to 
the same. The inner radiative corrections are included (see the text), but the outer radiative 
corrections are not. The long-dashed and dot-dashed lines are nearly indistinguishable in 
the lower panel. 



Inverse Beta Decay: Positron-Neutron Separation Vector 

The results above were based on lepton angular distributions as observed in 
water- Cerenkov detectors (light or heavy water). Some of the same reactions 
occur in scintillator detectors, but the lepton angular distributions are not 
observable due to the isotropic character of scintillation light. 

Despite that, thereJs a way to get directional information with a scintilla- 
tor detector. The ideaO makes use of the fact that scintillator detectors can 
measure final particle positions by the relative timing between different pho- 
totubes. For i/g + p — > e"*" -|- n, the positron is detected nearly at the point of 
creation. The neutron, however, will be detected (by its capture gamma rays), 
on average a few centimeters forward of the point of creation. The initially 
forward motion of the neutron is just a consequence of kinematics. Thus the 
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Figure 5: The average lepton cosine for the charged-current deuteron reactions, versus the 
neutrino or antineutrino energy, using Kubodera's calculations. We plot only from 1 MeV 
above threshold, so the De-dependence at low energies is not shown. The small jitters are 
due to the coarse integration grid. 



positron-neutron separation vector points in the incoming antineutrino direc- 
tion. 

However, before the neutron can be captured, it must be thermahzed. On 
the first several scatterings (about half of the kinetic energy is lost on each 
step), the forward motion tends to be preserved. In the remaining scatterings 
until thermal energy is reached, the scattering is isotropic. In the last phase, 
the neutron wanders randomly, undergoing further scatterings at thermal en- 
ergy, until it is captured. The neutron wander has the effect of degrading the 
significance of the initially forward motion. The average displ^icement and its 
fluctuations can be calculated by Monte Carlo simulationtilo. 

At the energies of reactor antineutrinos, the positron- neutron separation 
is about 1.5 cm and the statistical fluctiiation several cm, depending on the 
neutron capture time in the scintillator lj. There is an additional uncertainty 
in the position due to the gamma-ray localization error. For reactor antineutri- 
nos, the source direction is known, so a measurement of the average positron- 



18 



neutron separation vector can be used to make a background measurement, 
since background events will degrade the significance of the separation vector 
from the expectation. The precision with which the reactor can be located 
was evaluated by the Chooz collaboration and found to bea about 18°. Given 
the small displacement and the large error in localization, this is an impressive 
measurement. 

At the higher energies of supernova neutrinos, the neutron kinetic energy 
is larger and the elastic scattering cross section smaller, so the displacement 
is larger. The Chooz collaboration estimates^ that if they were to scale their 
detector to the size of SK (32 kton), that they would be able to locate a su- 
pernova to within about 9°, not so much worse than the neutrino-electron 
scattering result of about 5° in SK. The positron-neutron displacement tech- 
nique is not possible in SK since neutrons are not detected, and not possible 
in SNO because the neutron wanders over several meters. 

2.3 Triangulation 

Finally, there is the technique of triangulation by arrival-time differences. The 
idea is very simple - for two detectors, the cosine of the angle 9 between the axis 
connecting the two detectors (with separation d) and the supernova direction 
is determined by the arrival-time difference At by 

cos 61 = —. (22) 

For detectors on opposite sides of the Earth, d~AQ ms, and for SK and SNO, 
d = 30 ms. Two detectors thus define a cone of allowed directions to the 
supernova. Since the arrival time difference has an error (5(Ai), the cone will 
have an angular thickness 

(5(cos0) = ^^. (23) 

The timing error is obviously the crucial point, and in order to have 5[cos6) — 
0.1 (comparable to the neutrino-electron scattering result in SNO), one would 
need 5{/S.t) = 3 ms for the SK-SNO difference. This is much smaller than the 
duration of the supernova pulse (10 s) or even the risetime (assumed 100 ms). 

In order to assess the best that triangulation could do, consider the sce- 
nario where SK and the light water portion of SNO are each perfect detectors, 
identical except for size. (The shape of the event rate in the heavy water in 
SNO will be different, due to the different reactions, and it is not clear how 
to correct for that). Then the scattering rates observed in SK and the light 
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water in SNO will have the same underlying distribution, which depends on 
unknown details of the supernova models. 

The observed scattering rates in SK and the light water of SNO can then 
differ only in normalization, statistical fluctuations, and a possible delay. The 
normalizations matter only in that they determine the scale of the statistical 
fluctuations. SK will have about 10'' events and the light water in SNO about 
400 events, so the statistical error will be dominated by the timing error in 
SNO. That is, the underlying distribution will be measured in SK and this 
template will be used in SNO to test for a delay. 

In order to make estimates of the timing error, we need the shape of the 
event rate. We assume a short rise (in the current supernova models, this is 
somewhere between and 200 ms, so we will very conservatively use 100 ms), 
followed by a long decay (10 s, as for SN1987A). 

For the normalized event rate, we take 

f{t, to) = ai X — exp 

n 

f{t,tf)) = ^2 X — exp 

where 

ai = — ^ — , a2 = — - — . (26) 

Tl +T2 Tl + T2 

Then f{t, to) is a normalized probability density function built out of two 
exponentials, and joined continuously at t = to- In what follows, we assume 
that this form of /(i,to) is known to be correct and that ti and T2 are known. 
The shape of the event rate is shown in Fig. 6. 

The event rate in SNO will consist of N events sampled from f{t, to), and 
the event rate in SK will consist of N' events sampled from f{t,t'Q). Then 
At = to — t'o, and 6{At) ~ 6to, since the SNO error dominates. We consider 
only the statistical error determined by the number of counts. As noted, we 
want to determine the minimal error on the triangulation. This model, while 
simple, contains the essential timescales and an adjustable offset. 

These considerations lead to a well-posed statistical problem: If N events 
are sampled from a fcnoiiin distribution f{t, to), how well can to be determined? 
The Rao-Crapier theorem Ej provides an answer to this question (with a pos- 
sible subtletyLj). This theorem allows one to calculate the minimum possible 
variance on the determination of a parameter (here to), by any technique what- 
soever. This minimum variance can be achieved when all of the data are used 
as "efficiently" as possible, which is frequently possible in practice. 



{t - to) 

Tl 

jt-to) ' 

T2 



t < to 
t>to. 



(24) 
(25) 
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Figure 6: The schematic form of the normalized event rate f{t, to) is shown. To the left of to 
there is an exponential rise with time constant n . To the right of to there is an exponential 
decay with time constant T2 . The tick marks on the t-axis are in units of the respective time 
constants ri and T2. For clarity of display, we used ti/t2 — 10~^ in the figure, instead of 
the ri/r2 ~ 10~^ assumed in the analysis. 



One requirement of the theorem is that the domain of positive probability 
must be independent of the parameter to be determined. This condition is 
obviously not met for a zero risetime, since then the domain is (to, oo). For a 
nonzero risetime, the domain is technically (—00,00), independent of to, and 
so the theorem applies. The minimum possible variance on the determination 
of to is: 



1 



N X / dtf{t,to) 



din fit, to) 



dto 



(27) 



This is the general form for an arbitrary parameter to- When to is a translation 
parameter, i.e., f{t,to) depends only ont — to, this reduces to 



1 



Nx j dtf{t,to) 



din fit, to) 
dt 



(28) 
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= N X 



dt 



[df{t,to)/di 
f{tM) 



(29) 



For the particular choice of f{t,to) above, this reduces to 



1 



N X {ai/rl + aa/rl) , 



(30) 



('5^0 )L 



and so the minimum error is 



(-5^0) 



min 




(31) 



Note that Ni = N{ti/t2) is approximately the number of events in the ris- 
ing part of the pulse. Since the rise is the sharpest feature in f{t,to), it is 
unsurprising that it contains almost all of the information about tg. The to- 
tal number of events TV is fixed by the supernova binding energy release, so 
a change in the assumed total duration of the pulse, i.e., T2, would affect the 
peak event rate and hence the fraction of events in the leading edge. For a more 
general f{t,to), one would replace ti/t2 by this fraction computed directly. 

For SNO, iVi ~ 10-2 x 400 ~ 4, so 6{to) ~ 30 ms/^/4 ~ 15 ms. Since SK 
has about 25 times more events, the corresponding error would be about 3 ms. 
Therefore, the error on the delay is 5{At) ~ 15 ms and ^(cos6') ~ 0.50 at one 
sigma. We, have not specified the method for extracting to and hence At from 
the dataES. That is exactly the point of the Rao-Cramer theorem - that one 
can determine the minimum possible error without having to try all possible 
methods. 

IL-the risetime were zero, which seems to be unrealistic, then one can 



Here T2/N is simply the spacing between events near the peak. For a more 
general /{t^to), but still with a sharp edge at tg, one would simply replace T2 
by l/f(to). In fact, the shape of f{t,to) is irrelevant except for its effect on 
the peak rate, i.e., /(io)- So long as f{t,to) has a sharp edge and the right 
total duration, allowing a more general time dependence would therefore not 
change the results significantly. 

For the two cases, zero and nonzero risetime, we used different mathemat- 
ical techniques. This may seem like an artificial distinction, and that these two 
cases do not naturally limit to each other. In particular, it may seem incom- 
patible that in the first case the error ~ while in the second the error 




(32) 
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~ Further, one obviously cannot_take ri in the first result to obtain 

the second. However, it can be shownEa that the two techniques have disjoint 
regions of applicability (as a function of ri ) and that the formulas match nu- 
merically at the boundary between the two. The boundary is the value of ri 
such that the number of events in the rise is of order 1, i.e., the edge appears 
sharp for this or smaller ti. 

The final results for the pointing errors are given in Table IV. 

3 Concluding Remarks 

The next Galactic supernova should be a bonanza for neutrino physics and 
astrophysics. Should, that is, provided that we can make sense of the data. 
The difficulty is that there are many unknown aspects of both the particle 
physics properties of neutrinos and the astrophysics of the expected signal. 
Since tha. expected core-collapse supernova rate in the Galaxy is about 3 per 
centuryliJ, and since there are no neutrino detectors sensitive to supernovae in 
distant galaxies, we will not have the luxury of many observations. Of course, 
there is a worldwide effort to improve the numerical models of core-collapse 
supernovae l3. And besides neutrinos, there is one other direct probe of the 
supernova core, and that is gravitational radiation. LIGO may be ^nsitive to 
supernovae out to the Virgo Cluster (at rates of order 1 per yearE3), and its 
observations may improve the understanding of stellar collapse. 

Since there are uncertainties in the supernova models, the (t) test for the 
i^^ and Vt masses was designed to be as model-independent as possible. The 
resulting sensitivities for either j/^ or i^t are in Table III. If the and i^r are 
maximally mixed with a small mass difference, then the mass test is really on 
1712 = m^, and the limits in Table III improve by a factor of about %/2. 



Table 3; Mass sensitivity for various detectors and neutral-current reactions. The mass 
quoted is the limit which would be placed on either the or v-r mass at 90% CL if no delay 
were seen (equivalently, it is the smallest mass which would cause an unambiguous delaY)^ 
If the and Ut masses are nearly equal, then the sensitivity is improved by \/2. See also til 
for the -"^-^C results. 



detector 


neutral-current signal 


mass limit 


SNO 


u + d 


30 eV 


SK 





45 eV 


SK 


v + 


50 eV 


Kamland 




55 eV 


Borexino 




75 eV 
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There are also many possible signals of supernova neutrino oscillationsES. 
Ideally, the current and near-term terrestrial neutrino oscillation experiments 
can answer some of the key questions before the next supernova. Particularly 
crucial is the cojifistion of whether there are sterile neutrinos (or several flavors 
thereof). SNOEI may tell us whether the solar j/g neutrinos are oscillating into 
sterile or active flavors (or at all). Similarly for SK0, K2K0, and MINOSH 
and tha.atmospheric i/^ neutrinos. And MiniBoonellj may decide whether the 
LSNDcJ anomaly is correct (and hence whether sterile neutrinos are necessariil. 
These questions and more are discussed thoroughly in several recent reviewso. 

Since the neutrinos leave the supernova core a few hours before the light 
leaves the envelope, it should be possible for the neutrino experiments to give 
astronomers an advance warning of a Galactic core-collapse supernova, and also 
to provide some guidance as to where in the sky to look. Neutrino-electron 
scattering in SK seems to be the best technique, with an error of order 5°, de- 
pending on distance and the suppression of backgrounds from other reactions. 
While 5° may sound large to an astronomer, it is about 0.1% of the sky, which 
will make subsequent searches with small telescopes |iauch easier (the role of 
amateur astronomers is discussed in a recent article e3). Supernova location 
by triangulation seems to be rather difficult at present (see Table IV). 



Table 4: One-sigma errors on how well the direction to the supernova is defined by various 
techniques, at D = 10 kpc. The other parameters used are noted in the text. For neutrino- 
electron scattering, the most pessimistic background assumptions were used. 



technique 


error 


v + e'^ forward scattering (SK) 
V + e'^ forward scattering (SNO) 


59 ~ 5°, (5(cos6i) ~ 4 X 10"^ 
56 ~ 20°, (5(cos6i) ~ 6 x lO^^ 


j^e + P angular distribution (SK) 

Vp+p angular distribution (SNO) 

Ve d^Ve + d angular distributions (SNO) 


5{cos0) ~ 0.2 
5{cose) ~ 1.0 
5{cose) ~ 0.5 


triangulation (SK and SNO) 


5{cose) ~ 0.5 
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